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ABSTRACT 

Self-consistent treatment of time-dependent ionization in hydrodynamical calcu- 
lations of the X-ray emission from young supernova remnants has been performed. 
The novel feature of the calculations is that K a lines from species produced by 
inner-shell collisional ionization are included. Parameters of the shocked ejecta 
are found from fitting the model spectrum to the observed one. The application 
of the method to Tycho SNR using the classical deflagration model W7 for the 
explosion enables us to well reproduce the observed X-ray spectra and radial 
brightness profiles of the remnant. 

Key words: supernovae remnants - x-ray emission: time-dependent ionization 
- inner-shell ionization. 



1 INTRODUCTION 

The supernova remnant 1572 (Tycho SNR) is a con- 
venient object for testing thermonucle ar supernova 
mode ls. The availa bl e obs er vational data dSmi th et al."1 
| l98Sl: iFink et al. I Il994l : iDecourchelle et al. I l200lF 
IWem^net^hTTl27)05l) makes it possible to test both the 
theory of thermonuclear supernova explosions and gas- 
kinetic processes in pl asma heated by the SN blast wave. 
Stud i es of Tycho SNR dltoh et al. ll98Sl:lBrinkman et al. 1 
ll998l : ISorokina et al. ll2004l : iBadenes et al. 1120031) showed 
that it is difficult to find a model which describes 
self-consistently its origin and evolution. Nevertheless, 
the increasingly accurate account for various physical 
processes in analyzing of SNR observations makes it 
possible to narrow the range of possible explosion 

mechanisms. 

\ was shown in detail in ISorokina et al. I i2004t) 
(hereafter Paper I) , the accurate modeling of young SNRs 
even in one-dimensional hydro calculations requires tak- 
ing into account radiative energy losses and self-consistent 
calculation of plasma ionization state in every mesh zone 
at every time step. In addition, such processes as ther- 
mal conductivity of electrons and non-Coulomb energy 
exchange between electrons and ions should be included. 
It is very hard to treat these processes from the first prin- 
ciples, so two free parameters have been introduced. The 
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values of these parameters can be found from compari- 
son with observations. In Paper I we compared modeled 
emission with X-ray observations of Tycho SNR obtained 
by the XMM-Newton X-r ay observatory and reported in 
IDecourchelle et al. I (1200 lD . The data include X-ray (0.2- 
10 KeV) spectrum and azimuthally averaged profiles of 
X-ray brightness in silicon and iron lines (Fe XVII, Si K, 
Fe K). 

In Paper I numerical simulations of SN ejecta propa- 
gated into the surrounding interstellar medium have been 
performed for two models: the parametric deflagration 
model W7 llNomoto et al. Ill984l) . and new 3D (the "first 
principle" ) mod el mrO elaborated in M ax-Plank-Institute 
fur Astrophysik iReinecke et al. 120021) . X-ray spectra and 
radial brightness profiles of the remnant expected at an 
age of 430 years (the age of Tycho SNR) were computed 
for these models. However, the comparison with the XMM 
observations revealed noticeable discrepancies. For exam- 
ple, in the W7 model, X-ray spectrum did not reveal the 
presence Fe K line and for no model the observed relative 
position of X-ray brightness profiles in silicon and iron 
lines was obtained. However, these disagreements could 
not be used to reject the models considered, since not all 
relevant physical processes were fully incorporated into 
the radiation hydro co de at that time. 

In later studies bv lKosenko et af~l i2005l) , an attempt 
to modify the initial density profile of the SN ejecta was 
done. The outermost layers of the ejecta were made more 
rarefied. These modifications showed that the density pro- 
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file does not have a noticeable effect on the resulting X-ray 
brightness profiles of the remnant. However, a fine tun- 
ing of the parameters in the density distribution in the 
ejecta made it possible to well reproduce the observed ra- 
dial brightness profiles. In that case most of the emission 
in silicon and iron lines should have been generated by 
the external SNR shoc k (not by the reverse shock) , which 
seems rather unlikely llHwang et al. "l ll99Sl) . 

When calculating X-ray emission from the remnant 
in Paper I, we did not take into account inner-shell ion- 
ization processes from low-ionized elements. In station- 
ary plasma, collisional production of K a emission lines is 
suppressed due to a small value of the cross-section of in- 
teraction between free and inner-shell electrons. Usually 
K a emission lines from low-ionization ions are produced 
via photoionization processes (fluorescence). In shocked 
plasma time-dependent ionization occurs and in a SNR 
the ionization time-sc ale is about several hundred years 
dSorokina et al. 112004) . Thus the matter behind the re- 
verse shock front should be highly under-ionized com- 
pared to the temperature of ionizing electrons. Hence the 
collisional ionization of ions from inner-shells is highly 
operational i n the emission line production in a SNR. In 
IVink et al. I i2003Tl the role of inner-shell excitation pro- 
cesses to Si K a feature was investigated in SN1006. 

The main goal of this paper is to assess the role of 
inner-shell ionization processes in the X-ray emission of 
young supernova remnants. We calculate the hydrody- 
namic evolution and X-ray emission of a young remnant 
using the code SUPREMNA, introduced in Paper I, with 
inner-shell ionization processes taken into account. The 
example of the classical deflagration model W7 shows that 
it may produce the X-ray spectrum and radial brightness 
profiles very close to the ones observed from Tycho SNR 
for some set of parameters. In particular, the iron K a 
line may be pronounced in this model without invoking 
additional mixing in the ejecta. 



2 CALCULATION OF THE INNER-SHELL 
IONIZATION IN SNRS 

To calculate the intensity of a line produced by ioniza- 
tion from the inner shell (n = 1) of an ion followed by 
the transition of an electron from level n to this vacant 
place, the following fo rmula was employed JSpitzer Il978l: 
iKallman et al. 112004) 



Jnl = Hi n e {<Jv)i UJi. n l E n l 



r erg 



(1) 



lcm° sJ 

where rii, n e are number densities of ions i and electrons 
correspondingly ; 

{' oo 

(av)i = / avf(v)dv (2) 

■>E tht [i) 

is the velocity averaged cross-section. Here -Ethr[i] is the 
threshold ionization energy from the inner-shell of ion 
i, f(v) is the Maxwellian distribution function for elec- 
trons; uii.ni is the fluorescence yield (photons/ionization) 
for the line which corrects the intensity for photonless 
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Figure 1. Ionization rates for Ar,Ca,Fe,Ni. Asterisks ("*") 
show standard ionization rates used in SNR simulations. Cir- 
cles ("o") show the inner-shell ionization rates. Matter with 
an initial temperature of T; n = 10 4 K was instantly heated 
up to Tfj n = 2 X 10 s K. The ionization parameter is nt = 
3 X 10 9 s/cm 3 . 



transitions due to Au ger electrons llKaastra. Mewe Il993t 
IKallman et al. ll2004) : E n i is the energy of the line. 

Approximate formulas for the inner-shell io n izatio n 
cross-sections a were taken from iHombourger I lll99ctl . 
These formulas fit experimental data for a wide range of 
ions and electron energies. Each cross-section formula was 
convolved numerically with the Maxwellian velocity dis- 
tribution (formula[5Jl. For each species dependence (av)i 
on temperature was numerically approximated with an 
accuracy of better than 1% and then inserted into the 
emis sion routine. Da t a on E t hi [i] , Wj.ni , E n i was taken 
from iKaastra. Mewe I (ll99St) . 

In the SNR modeling, the plasma ionization state is 
calculated by solving a system of kinetic equations in each 
grid mesh at every time step. Assuming that inner-shell 
ionization processes do not affect very much the electron 
production rate and thus the equation of state of the mat- 
ter, the inner-shell ionization has not been included in the 
kinetic equations. Fig0 shows ionization rates for some 
ions. Ionization rates which are included in the present 
paper (and in Paper I ) are shown by asterisks "*" , the 
inner-shell ionization rates are presented by circles "o". 
In most cases the inner-shell ionization rates are several 
times smaller than the standard (outer-shell) ionization 
rates. 

We did not include inner-shell excitation processes 
in our method. Contribution of excitations from inner- 
shell is negligible compared to normal excitations and 
inner-shell ionization rates. For example, formulae from 
iMewe. Gronenschild I (|l98ll ) show that the ratio of nor- 
mal excitation rates to inner-shell ionization and to inner- 
shell excitation for Fe XXV at temperature ~ 10 s K is 
about 1 : 0.5 : 10" 4 . 
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3 COMPARISON WITH OBSERVATIONS 

The results of the calculations have been compared with 
the X-ray observations of Tycho SNR performed by 
XMM-Newton observatory. The data were taken from the 
public available library of the observatory 1 . These obser- 
vations were conducted with the EPIC instrument, the 
spectrum was extracted from data obtained by the MOS1 
camera. 

Based on W7 explosion model, eight variants of hy- 
drodynamical calculations have been performed for differ- 
ent values of three parameters describing the thermal con- 
ductivity (Ckiii), non-Coulomb energy exchange between 
electrons and ions (q) and the density of the circumstellar 
matter (pcsm). Calcul ated spectra were converted to an 
XSPEC llArnaud |ll99rf) table model with these three in- 
terpolation parameters. The physical meaning of the pa- 
rameters Ckm and q is as follows. 

The factor < Ckm < 1 describes the fraction of 
the thermal flux contributed to the system compared to 
the standard case (formula (7) in Paper I), where the 
flux is limited only by the existence of a maximum speed 
(the speed of sound) for the heat carriers, ignoring a 
possible decrease in the particle mean free path due to 
magnetic fields and plasma instabilities. To describe the 
effects of collisionless energy exchange, in Paper I we 
introduced the parameter q that specifies the fraction 
of artificial viscosity Q, added to the pressure of ions: 
P% = Pi (thermal) + qQ. Then for the electronic compo- 
nent we put P e = P e (thermal) + (1 — q) Q. If only the col- 
lisional exchange is taken into account, q = 1 and we used 
the standard system of equations with only ions heating 
of at the front. 

Comparison of the spectra 

By fitting the XSPEC table models to the XMM- 
Newton observations, we found the values of the table 
interpolation parameters. Then with these values a new 
hydrodynamical model was calculated. Fig|2] shows the 
observed spectrum (EPIC MOS1) and theoretical spec- 
trum obtained for the W7 model. Table Q lists the values 
of fitting parameters. The best-fitting values for the col- 
umn density Nh is also shown in the table. 

Inset in Fig|U shows Fe K line profile in details. It 
is clear from the figure that modeled centroid shifted to 
higher energies, compared to observable one. This discrep- 
ance indicates that the matter in our model is slightly 
overionized, thus the temperature or ionization timescale 
of the Tycho remnant should be lower. These arguments 
point to less energetic model for the remnant. 

Comparison of the remnant's brightness pro- 
files 

Following the scheme elaborated in Paper I, the ra- 
dial brightness profiles of the model remnant were con- 
struc ted in Si K, Fe XVII, Fe K lines .Decourchelle et al~l 
l200lT) . The results are plotted in Fig|3] The figure shows 
the X-ray brightness profiles for W7 model, calculated 



1 http://xmm.vilspa.esa.es/external/xmm_data_acc/xsa/ 
obsID = 009621010 
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Figure 2. Observed (the gray crosses) and theoretical (the 
black solid line) spectra of the Tycho SNR. Inset shows Fe K a 
line in details. 



Table 1. The values of the best fit parameters for Tycho SNR 
model 



PCSM, g/cm 3 q Ckm N H , cm 



1.8 x 10" 24 0.93 0.0085 3.7 x 10 



with the best-fitting interpolation parameters (Table 0. 
These profiles are very similar to the ones found in 
iDecourchelle et al. I l)200ll) . 

Using the model radial profile, the distance to the 
remnant was evaluated. Table [5] shows estimates of the 
distance to Tycho SNR as obtained from the X-ray flux 
normalization (column 2) and as derived from the bright- 
ness profiles (column 3). Also the value of distance to 
Tycho remna t from observed proper motion of filaments 
ISmith et al. I dl99lT) presented. The distance estimated 
from normalizations of the SNR spectrum is more than 
two times smaller than that derived from the brightness 
profiles. This may suggests that our model of Tycho SNR 
is underluminous. 

Table 2. The estimates of the distance to Tyc ho SNR from 
W7 m odel, and estmimates of the distance from ISmith et all 
Il99lf> 



flux 


profiles 


H a filaments 


1.34 ±0.01 kpc 


3.1 ±0.1 kpc 


1.5 - 3.1 kpc 
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Figure 3. Simulated brightness profiles for the W7 model. 
The solid black line shows the density profile. Step-like profiles 
are observed by XMM-Newton, taken from Decourchelle et al. 
(2001): Fe K (solid line), Fc XVII (dashed line). 



4 DISCUSSION 

In this paper we have performed a hydrodynamical cal- 
culation of X-ray spectrum and radial brightness profiles 
in X-ray lines of SNR 1572 and compared the model X- 
ray emission with the XMM-Newton observations. We 
have used the classical deflagr ation model of therm onu- 
clear supernova explosion W7 dNomoto et al. lfl984l) . We 
employed the method presented in Paper I and comple- 
mented by the inner-shell collisional ionization processes 
in the calculations of X-ray emission from the remnant. 
This improvement allowed us to remove discrepancies be- 
tween theoretical and observational data found in Paper 
I. 

We compared the model spectrum with the observ- 
able one by constructing the xspec table model with 3 
interpolation parameters: q, Ckm, Pcsm- This model was 
used to fit observational data and to find the best-fitting 
values of the interpolation parameters. The moderate 
value of the parameter Ckm — 0.01 and moderate value 
of the parameter q ~ 0.9 suggest that thermal conduc- 
tion is not very important and points to an appreciable 
non-Coulomb energy exchange between electron and ion 
components. These values also suggest the presence of 
magnetic field and plasma instabilities in the ejecta of 
the SNR. 

The derived column densit y Nh is slightly smalle r 
than that es t imated earlier b y |Albinsoji_et L jaLj tl986|). 
ISmith et all lll98Sl) . Iltoh et al. I (Il988t ). iHwang et al. 1 



Jl998f) an d the value of pcsm is hi gh er than t hat de - 
rived bv iDecourchelle et al. I feOOll) . iHuehes I J2000|) . 
lBrmkmanet^u^ " j^98ti '. Since it is not excluded that the 
W7 model is not the best one for the Tycho SNR, these 
parameters should be considered as zero approximation. 
In forthcoming studies we plan to examine other mod- 
els with different modes of burning in type la supernovae 
with the goal to select the best-fitting one. 

Estimations of the distance to the remnant using 
the X-ray flux normalization and remnant's angular size 
(Table disagree likely due to neglecting the nonther- 
mal component from the forward shock in X-ray emission 
modeling. Furthermore, probably the explosion model of 
Tycho supernova was more energetic than W7 model. 



5 CONCLUSIONS 

Our conclusions are as follows. The calculated X-ray spec- 
trum and radial brightness profiles indicate that even the 
classical SN la deflagration model W7 is capable to re- 
produce the main observable features of a young SN la 
remnant, provided that all important physical processes 
are taken into account. In this case there is no need to 
introd uce additional elem ent mixing in the ejecta of this 
model (lltoh et al. II198&J) 
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